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Abstract 

Recently, the shell model in the complex /c-plane (the so-called Gamow Shell 
Model) has been formulated using a complex Berggren ensemble representing 
bound single-particle states, single-particle resonances, and non-resonant con- 
tinuum states. In this framework, we shall discuss binding energies and energy 
spectra of neutron-rich helium and lithium isotopes. The single-particle basis 
used is that of the Hartree-Fock potential generated self-consistently by the 
finite-range residual interaction. 
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I. INTRODUCTION 



Low-energy nuclear physics is undergoing a revival with revolutionary progress in ra- 
dioactive beam experimentation. New facilities have been built or are under construction, 
and new ambitious future projects, such as the Rare Isotope Accelerator in the U.S.A., 
will shape research in this field for decades to come. From a theoretical point of view, 
the major problem is to achieve a consistent picture of weakly bound and unbound nuclei, 
which requires an accurate description of the particle continuum properties when carrying 
out multi-configuration mixing. This is the domain of the continuum shell model [1] and, 
most recently, the Gamow Shell Model (GSM) [2,3] (see also [4]). GSM is the multicon- 
figurational shell model with a single-particle (s.p.) basis given by the Berggren ensemble 
[5] which contains Gamow (or resonant) states and the complex non-resonant continuum. 
The resonant states are the generalized eigenstates of the time-independent Schrodinger 
equation which are regular at the origin and satisfy purely outgoing boundary conditions. 
The s.p. Berggren basis is generated by a finite-depth potential, and the many-body states 
are obtained in shell-model calculations as the linear combination of Slater determinants 
spanned by bound, resonant, and non-resonant s.p. basis states. Hence, both continuum 
effects and correlations between nucleons are taken into account simultaneously. An inter- 
ested reader can find all details of the formalism in Refs. [2,3] in which the GSM was applied 
to many-neutron configurations in neutron-rich helium and oxygen isotopes. In this con- 
tribution, we shall present the first application of the GSM formalism to the p-shell nuclei 
in the model space involving both neutron and proton s.p. Gamow states calculated from 
the self-consistent Hartree-Fock (HF) method (Gamow Hartree-Fock, GHF). Moreover, we 
briefly report on the first successful application of the density matrix renormalization group 
(DMRG) technique [6] in the context of realistic multiconfigurational shell model (SM). 



In our previous studies [2,3], we have used the Surface Delta Interaction (SDI) and the 
s.p. basis has been generated by a Woods-Saxon (WS) potential which was adjusted to 
reproduce the s.p. energies in 5 He. This potential (" 5 He" parameter set [3]) has the radius 
R = 2 fm, the diffuseness d =0.65 fm, the strength of the central field V = 47 MeV, and the 
spin-orbit strength V so =7.5 MeV. The SDI interaction has several disadvantages in practical 
applications. Firstly, it has zero range, so one is forced to introduce an energy cutoff and, 
consequently, the residual interaction depends explicitly on the model space. Moreover, one 
is bound to use the same WS basis for all nuclei, as the SDI interaction cannot practically be 
used to generate the HF potential. Consequently, as the chosen WS basis is not an optimal 
s.p. basis (HF basis is), one cannot easily truncate the configuration space when the number 
of valence particles increases [3] . So, we have decided to introduce a new two-body residual 
interaction, the Surface Gaussian Interaction (SGI): 



II. DESCRIPTION OF THE CALCULATION 
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which is used together with the WS potential with the " 5 He" parameter set. 
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The SGI interaction is a compromise between the SDI and the Gaussian interaction. The 
parameter R in Eq. (2.1) is the radius of the WS potential, and Vq(J,T) is the coupling 
constant which explicitly depends on the total angular momentum J and and the total 
isospin T of the pair of nucleons. A principal advantage of the SGI is that it is finite-range, so 
no energy cutoff is needed. Moreover, the surface delta term in (2.1) simplifies the calculation 
of two-body matrix elements, because the radial integrals become one-dimensional and they 
extend from r=0 to r=2R$. (In the Gaussian case, such as the Gogny force, they are 
two-dimensional and have to be extended to infinity.) Consequently, an adjustment of the 
Hamiltonian parameters becomes feasible. Finally, the resulting spherical HF potential 
is continuous and can be calculated very accurately. This allows one to use the optimal 
spherical HF potential for the generation of the Berggren basis for each nucleus studied; 
hence a more efficient truncation in the space of configurations with a different number of 
particles in the non-resonant continuum. In the present study, it turned out to be sufficient 
to consider at most two particles in the GHF continuum. This restriction on the number 
of particles in the non-resonant continuum allowed us to extend the studies up to unbound 
10 He and the halo nucleus 11 Li. 



A. Choice of the valence space 

In our He and Li calculations, the valence space for protons and neutrons consists of 0^3/2 
and Opi/2 spherical GHF resonant states, calculated for each nucleus, and the {^3/2} and 
{ipi/2} (i — 1, • • • , M) complex continua generated by the same potential. These continua 
extend from i?e[A;]=0 to i?e[fc]=8fm _1 and they are discretized with 14 points (i.e., M = 14). 
Altogether, we have 15 JO3/2 and 15 pi/ 2 GHF states (shells) in the GSM calculation. The 
imaginary parts of k- values of the discretized continua are chosen to minimize the error made 
in calculating the imaginary parts of energies of the many-body states. Other continua, such 
as S1/2, d 5 / 2 , ■ ■ ■ are neglected, as they can be chosen to be real and would only induce a renor- 
malization of the two-body interaction. We have checked [2,3] that their influence on the 
binding energy of light helium isotopes is negligible. On the other hand, the lsi/2 anti-bound 
neutron s.p. state is important in the heaviest Li isotopes ( 10 Li, n Li) and plays a significant 
role in explaining the halo ground-state (g.s.) configuration of n Li [7,8]. At present, however, 
solving a GSM problem for n Li in [0p 3/2 , {^3/2}; OP1/2, {ipi/2}] lsi/2, {^1/2}; (i = 1, • • • , M)\ 
GHF space is not possible within a reasonable computing time. This task will be, how- 
ever, possible in the near future by using a new generation GSM code which employs the 
DMRG methods to include the non-resonant continuum configurations contribution in the 
many-body wave function [9,10]. 

Having defined a discretized GHF basis, one constructs the Slater determinants from all 
s.p. basis states (bound, resonant, and non-resonant), keeping only those with at most two 
particles in the non-resonant continuum. Indeed, as the two-body Hamiltonian is diagonal- 
ized in its optimal GHF basis, the weight of configurations involving more than two particles 
in the continuum is usually quite small, and they are neglected in the following. 
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B. The Helium chain 



Within the chain of helium isotopes, which are described assuming an inert 4 He core, 
there are only T—l two-body matrix elements. Consequently, only (J=0, T—l) and (J=2, 
T—l) couplings come into play. We have adjusted V (J — 0, T — 1) to reproduce the 
experimental g.s. energy of 6 He relative to the g.s. of 4 He, whereas Vq(J — 2,T — 1) has 
been fitted to all g.s. energies from 7 He to 10 He. Indeed, these latter states are mainly 
sensitive to V$(J = 2,T = 1) whereas, for obvious geometrical reasons, the J = 2,T = 1 
coupling is absent in the g.s. of 6 He. The adopted values are: V (J — 0, T — 1) = -403 
MeV fm 3 and V (J = 2, T = 1) = -315 MeV fm 3 . 

The calculated spectrum of He isotopes is shown in Fig. 1. The experimental g.s. binding 
energies relative to the 4 He core are reproduced fairly well. For instance, the g.s. of 6 Hc 
and 8 He are bound whereas 5 He and 7 He are unbound. Moreover, the so-called 'helium 
anomaly' (see, e.g., Ref. [15]), [i.e., the higher one- and two- neutron emission thresholds 
in 8 He than in 6 He], is well reproduced. The qualitative features seen in the experimental 
spectra are also satisfactorily reproduced by the GSM with the SGI Hamiltonian. The width 
of certain resonance states is too large mainly because these states are calculated at too high 
an excitation energy above the one- and two-nucleon emission threshold. 

Looking at the configuration mixing in the g.s. of helium isotopes (see Tables I-V), one 
can see that the coupling to the non-resonant continuum is most important in 6 He and 7 He. 
This finding is consistent with the Borromean nature of 6 He, which is bound only because of 
correlation effects (pairing scattering to the continuum). The g.s. of 7 He is not a simple s.p. 
resonance, and this is reflected in GSM wave functions by large amplitudes of configurations 
involving particle(s) in the non-resonant continuum. This in turn implies that the g.s. width 
of 7 He is a result of the complicated mixture of different configurations with three valence 
neutrons in various resonant and non-resonant shells. 

The coupling to the continuum is somewhat less important in the heavier isotopes ( 8 He, 
9 He, and 10 He). This is partially related to the fact that these nuclei have a closed 
shell; hence for the bound g.s. of 8 He, the calculated spherical HF potential is the exact HF 
potential (no deformation effects are present). Consequently, the contributions of the lp- 
lh excitations to the g.s. wave function vanish. For unbound states, there are small lp-lh 
components due to the neglect of the imaginary part in the approximate HF potential used in 
the actual calculation. Secondly, in the heavier He isotopes, the GHF Ops/2 shell is calculated 
to be bound, which further diminishes the importance of the coupling to the non-resonant 
continuum. Nevertheless, this coupling is by no means negligible, as the configurations 
involving states of the non-resonant continuum still represent about 10% of the wave function 
and play a significant role in generating binding for all these nuclei. 

Comparing the SGI results with those of Ref. [3], one can see that the main conclusions 
for 6 He and 7 He still hold. This is because the HF potential for these nuclei is close to the WS 
potential which was previously used. However, once the 0p 3 / 2 shell gets closed, the HF field 
becomes very different from the WS potential. As a consequence, the large configuration 
mixing obtained in Ref. [3] in the WS basis for 8 He is significantly reduced when the GHF 
basis is used. In other words, a large part of the configuration mixing in the WS basis is 
not due to the presence of genuine two-particle correlations involving continuum states, but 
due to the lp-lh couplings now incorporated in the GHF basis. 
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In order to assess the importance of two-particle correlations involving continuum states, 
we have calculated the g.s. of 6 He, 7 He, and 8 He, but in the model space of lp-lh excitations 
only. In such truncated calculations, the real part of the energy reads: 0.75 MeV for 5 He, 

0. 58 MeV for 6 He, 1.27 MeV for 7 He, and -1.40 MeV for 8 He. Based on these numbers, 
one is tempted to conclude that the anomalous increase of one- and two-neutron separation 
energies in 8 He, as compared to 6 He, is a genuine mean-field effect caused by the J = 2, T = 1 
coupling. The coupling to the non-resonant continuum further enhances the helium anomaly. 

The drip-line nucleus 8 He is two-neutron bound already in the HF approximation, so 
the coupling to the non-resonant continuum provides only additional binding. On the other 
hand, 6 He is unbound in GHF, so the continuum coupling is solely responsible for binding. 
Therefore, one may conclude that the Borromean features in the 4_8 He chain are caused 
mainly by the continuum couplings, whereas the helium anomaly is mainly due to the 
change from a ( J = 0, T = l)-dominated mean field to a mean field dominated by the 
(J = 2,T= 1) coupling. 

C. The Lithium chain 

As a first example of GSM calculations in the space of proton and neutron states, we have 
chosen to investigate the Li chain. The continuum effects are very important in these nuclei, 
both in their ground states and in excited states. The nucleus n Li is also a well-known 
example of a two- neutron halo. In our p-space(s) calculation, we consider the one-body 
Coulomb potential of the 4 He core, which is given by a uniformly charged sphere having 
the radius of the WS potential. It turns out that the inclusion of the one-body Coulomb 
potential modifies the GHF basis in lithium isotopes as compared to the helium isotopes, 
an effect which is usually neglected in the standard SM calculations. 

Recent studies of the binding energy systematics in the s<i-shell nuclei using the Shell 
Model Embedded in the Continuum (SMEC) have reported a significant reduction of the 
neutron-proton T=0 interaction with respect to the neutron-neutron T—l interaction [16,17] 
in the nuclei close to the neutron drip line. In SMEC [18], this reduction is associated with 
a decrease in the one- neutron emission threshold when approaching the neutron drip line, 

1. e., it is a genuine continuum coupling effect. The detailed studies in fluorine isotopes have 
shown that the reduction of the T=0 neutron-proton interaction cannot be corrected by 
any adjustment of the monopole components of the effective Hamiltonian. To account for 
this effect in the standard SM, one would need to introduce a iV-dependence of the T=0 
monopole terms. Interestingly, it has recently been suggested [19] that a linear reduction of 
T=0 two-body monopole terms is expected if one incorporates three-body interactions into 
the two-body framework of a standard SM. 

Our GSM studies of lithium isotopes indicate that the reduction of T=0 neutron-proton 
interaction with increasing neutron number is essential. For example, if one uses the 
Vo(J, T = 0) strength adjusted to 6 Li to calculate 7 Li, the g.s. of 7 Li becomes overbound 
by 13 MeV, and the situation becomes even worse for heavier Li isotopes. To reduce this 
disastrous tendency, in the first approximation we have used a linear dependence of T=0 
couplings on the number of valence neutrons n: 

V ( J = 1, T = 0) = a w [1 - /?i (n " 1)] . 
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V (J = 3, T = 0) = a 3 o [1 - fcoin - 1)] , 



(2.2) 



with a 10 = -600 MeV fm 3 , j3 w = -50 MeV fm 3 , a 30 = -625 MeV fm 3 , and (3 30 = -100 
MeV fm 3 . This linear dependence is probably oversimplified, as shown in Refs. [16,17] where 
the proton-neutron T=0 interaction first decreases fast with increasing neutron number and 
then saturates for weakly bound systems near the neutron drip line. For the T—l interaction, 
we have taken the parameters V (J — 0,T — 1) and V (J — 2,T — 1) determined for the 
He chain (see Fig. 1 and the preceding discussion). 

The results of our GSM calculations for the neutron-rich Li isotopes are shown in Fig. 2. 
One obtains a reasonable description of the g.s. energies of lithium isotopes relative to the 
g.s. energy of 4 He, even though the agreement with the data is somewhat worse than in the 
He chain. Clearly, the particle-number dependence of the T=0 matrix elements has to be 
further investigated. The absence of an antibound Si/ 2 state in the Berggren basis is most 
likely responsible for large deviations with the data seen in 10 Li and 11 Li. 

III. FUTURE PERSPECTIVES: THE DENSITY MATRIX RENORMALIZATION 
GROUP TECHNIQUES FOR SOLVING THE GAMOW SHELL MODEL 

As outlined above, in the GSM one uses a Berggren basis which consists of discrete 
states (bound and resonant states, I = anti-bound state for neutrons) and the discretized 
non-resonant continuum. Consequently, the dimension of the (non-hermitian) Hamiltonian 
matrix in GSM grows extremely fast with increasing the size of the Hilbert space. This 
'explosive' growth of the dimension is much more severe than in the standard SM for which 
the dimensionality problem concerns only the 'pole' space of GSM. The future perspectives 
of GSM applications are ultimately related to the progress in developing new methods of 
truncating huge SM spaces. One promising approach is the DMRG method. In nuclear 
structure, this method has been successfully applied in schematic Hamiltonians [21], but no 
fully convincing results have so far been reported in the context of the realistic SM. 

Let us begin by summarizing the basic elements of DMRG [6] (see also Ref. [21] for a 
pedagogical discussion relevant to nuclear physics applications). The main idea is to consider 
'step by step' different s.p. shells in the configuration space and retain only the N opt "best 
states" dictated by the one-body density matrix. The convergence of the DMRG method is 
then studied with respect to N opt . 

The configuration space is divided into two subspaces denoted by 7i and V . Generally 
speaking, H contains the lowest s.p. shells. In the first step, one calculates and stores all 
the possible matrix elements of suboperators of the Hamiltonian in the Ti space: 

a\ (a*a) j , (aV^, ((aW) j a) J ), ((aV) J '(5a) J ')). 

One also constructs in H all states \h) with 0, 1,2, • • • particles coupled to all possible J- 
values. Then, one considers the first s.p. shell in V, calculates matrix elements of all 
suboperators in this shell, and constructs in V all states \p) with 0, 1, 2, • • • particles coupled 
to all possible J-values. In the following, one adds 'one by one' additional s.p. shells in 
V: one calculates all matrix elements of suboperators in the added shell and all states with 
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0, 1, 2, • • • particles in the added shell and in the shell previously considered in V. New states 
are successively added in V until the number of states \p) with 0, 1, 2, • • • particles is larger 
than N opt . Then one diagonalizes the Hamiltonian in the space \hp) J made of vectors in 
H and V. Obviously, the number of particles in such states is equal to the total number 
of valence particles in the system, and J is equal to the angular momentum of the state of 
interest. From the eigenstates : 



The density matrix is diagonalized and only N opt eigenstates having the largest eigenvalues 
are retained. (In the non-hermitian GSM problem, eigenvalues of the density matrix are 
complex and the eigenstates are selected according to the largest absolute value of the density 
eigenvalues.) One then recalculates all the matrix elements of suboperators for the optimized 
states; they are linear combinations of previously calculated matrix elements. Then, one 
adds the next shell in V and, again, only the N opt states selected according to the above 
prescription are kept. This procedure is continued until the last shell in V is reached, 
providing a 'first guess' for the many-body wave function. 

At this stage, a "sweeping phase" begins in which the iterative process is reversed. First, 
one considers the last shell in V. At this point, one constructs states with 0, 1, 2, • • • particles 
and then the process continues, step by step, until the number of vectors becomes larger 
than N opt . When the i th shell in V is reached, the Hamiltonian is diagonalized in the set of 
vectors \h,p prev ,p) J , where h is a state that belongs to Tl, p prev is a previously optimized 
state first shells in V), and p is a state that concerns all the shells between the i th one 
and the last one in V . The density matrix is then diagonalized and the N opt p-states are 
kept. The procedure continues by adding the (i — l) st shell in V, etc., until the first state 
is reached. Then the procedure is reversed: the first shell is added, then the second, the 
third, etc. The succession of sweeps is successful if the energy gradually converges after 
every sweep. 

As a first application of DMRG, we calculated the g.s. 0| and the excited 0^~ state in 6 He 
in the same configuration space as in Sec. II B. Initially, we applied the DMRG procedure 
without sweeping, i.e., in its 'infinite algorithm' version [6]. In this calculation, we considered 
14 pi/2 shells and 14 p 3 / 2 shells in the non-resonant continuum. The method succeeded in 
reproducing the "exact" g.s. energy of 6 He when taking 16 shells in H (8 pi/ 2 shells and 8 
p 3 /2 shells) and keeping 6 vectors in V at each iteration. One should mention that the states 
in 7i have not been optimized during the iterative procedure. This example demonstrates 
that the DMRG in the infinite-system variant can be generalized for genuinely non-hermitian 
problems. On the other hand, the resulting gain in reducing the dimensionality of the GSM 
is not particularly impressive. 

The gain factor is radically improved when the finite-system algorithm (sweeping) is 
applied. In this case, for 20 pi/ 2 shells and 20 p 3 / 2 shells in the non-resonant continuum, 
an excellent convergence has been reached for both states of 6 He by taking 0pi/ 2 and 0p^/ 2 
Gamow resonances in TL and keeping in V only six vectors after adding each shell. In the 



(3.1) 



one calculates the one-body density matrix: 




(3.2) 
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considered example, the total dimension of the GSM Hamiltonian is 462 and the rank of the 
biggest matrix to be diagonalized in GSM+DMRG is 32. The gain factor is expected to be 
even more impressive for a larger number of valence particles. 

IV. CONCLUSIONS 

The Gamow Shell Model, which has been introduced only very recently [2,3], has proven 
to be a reliable tool for the microscopic description of weakly bound and unbound nuclear 
states. In He isotopes, GSM with either SDI or SGI interactions was able to describe 
fairly well the many-body properties, in particular the Borromean features in the chain 
4_8 He. Using the finite-range SGI interaction made it possible to perform GHF calculations, 
thus designing the optimal Berggren basis for each nucleus. In this way, we were able to 
disentangle the correlations due to the continuum coupling from the particle-hole excitations 
essential for building the mean field. We have shown that the Borromean features of the 
helium isotopes are the results of the correlations involving the non-resonant continuum, 
whereas the helium anomaly is a mean field effect due to the transition from a(J = 0,T=l)- 
dominated mean field in 6 He to the ( J = 2, T = 1) GHF field in 8 He. 

In Li isotopes, T=0 matrix elements of the two-body interaction could be studied for 
bound and resonant many-body states. It was found that the T=0 interaction contains a 
pronounced density (particle-number) dependence which originates from the coupling to the 
continuum and leads to an effective renormalization of the neutron-proton coupling. This 
effect cannot be absorbed by the modification of T=0 monopole terms in the standard SM 
framework. The effective renormalization of (J = 1,T = 0) and (J = 3, T — 0) couplings, 
which has been found in the present GSM studies, has to be further investigated, as it is also 
related to the question about the importance of 3-body correlations and density dependence. 

The successful application of GSM to heavier nuclei is ultimately related to the progress 
in optimization of the GSM basis. The promising development, discussed in this paper, is 
the adaptation of the DMRG method [6] to the genuinely non-hermitian SM problem in the 
complex-A; plane. The first applications using the j-scheme GSM are very promising; they 
demonstrate that we may well be at the edge of solving the GSM in today's inaccessible 
model spaces. 

This work was supported in part by the U.S. Department of Energy under Contract 
Nos. DE-FG02-96ER40963 (University of Tennessee) and DE-AC05-00OR22725 with UT- 
Battelle, LLC (Oak Ridge National Laboratory), DE-FG05-87ER40361 (Joint Institute for 
Heavy Ion Research), and by the Polish Committee for Scientific Research (KBN). 
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TABLE I. Squared amplitudes of configurations in the ground state of 6 He. The sum of squared 
amplitudes of Slater determinants with n particles in the continuum is denoted by . 



Configuration 


c 2 


Qp? /2 

r(l) 
L + 


0.656-i0.566 
6.06-10~ 3 -i0.0516 
0.363+10.509 


r(2) 
L + 


-0.0245+i0.108 


TABLE II. Similar as in Table I except for the 3/2]" g.s. of 7 He. 


Configuration 


c 2 


0^3/2 

0pJ /2 P 2 1/2 

0^3/2 PPl/2 

r(l) 
^+ 

r(2) 
^+ 


0.331-10.0973 
0.0111-i0.0406 
5.380-10- 4 -i6.950-10~ 4 
0.507+10.0714 
0.150+i0.0672 
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TABLE III. Similar as in Table I except for the g.s. 


of 8 He. 


Configuration 


c 2 


°Pt/2 

°P 2 3/2 PP?/2 

r(l) 
^+ 

r(2) 
^+ 


0.889-i7.826-10~ 3 
0.0316-10.0529 
0.0613+10.0226 
0.0184+10.0225 




TABLE IV. Similar as in Table I except for the 1 /2\ \ 


l.s. of 9 He. 


Configuration 


c 2 


0^3/2 °P{/2 

r(l) 
^+ 

r(2) 
^+ 


0.937-10.0209 
0.0473-il.779-10~ 3 
0.0157+10.0227 




TABLE V. Similar as in Table I except for the g.s. 


of 10 He. 


Configuration 


c 2 


0^3/2 °Pl/2 

r(l) 
^+ 

r(2) 
^+ 


0.965-10.0206 
-5.640-10- 3 +i8.392-10~ 3 
0.0409+i0.0122 
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FIG. 1. GSM spectra of helium isotopes obtained with the SGI Hamiltonian. Experimental 
data are taken from Refs. [11-14]. 
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FIG. 2. GSM spectra of lithium isotopes obtained with the SGI Hamiltonian. Experimental 
data are taken from Refs. [11,12,20]. 
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